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1.  Introduction 

This  paper  is  concerned  with  derivation  of  an  optimal  decision  rule 
from  an  integral  representing  the  costs  associated  with  adjustments  in 
production  to  changing  sales  forecasts.  The  cost  formulation  used  here 
and  the  optimal  rule  associated  with  it  are  closely  related  to  that 
examined  by  Holt,  Modigliani,  Simon,  and  the  present  writer  [2,1].  However, 
the  simple  case  to  be  examined  will  allow  some  of  the  dynamic  performances 
characteristics  of  the  scheduling  rules  and  the  effects  of  errors  in 
estimates  and  expectations  to  be  examined  in  a  relatively  convenient  manner. 


*In  the  preparation  of  this  paper  I  havo  benefited  from  discussions  with 
Charles  C.  Holt,  Franco  Modigliani,  ana  Herbert  A.  Simon. 
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As  a  first  approximation,  the  costs  of  inventory  storage  and  depletion 
and  of  overtime  pay  to  the  work  force  of  a  factory  might  be  represented  by 
the  following  integral: 

(1.1)  C  -  lim  2<h“h  )2  ♦  (p-P  )2  dtl 

T-o>  ^  -T  L  c  C  J 

where  C  is  the  cost  per  unit  of  time  associated  with  inventory  and  overtime 
decisions,  7 \  represents  the  relative  costs  of  decx*easing  inventories  to 
increasing  the  amount  of  overtime  of  the  woifc  force,  Hc  and  Pc  are  the 
minimum-cost  inventory  level  and  production  rater  respectively,  H(t) 
represents  the  level  of  finished- goods  inventories  at  time  t,  P(t)  represent? 
the  production  rate  at  time  t,  and  S(t)  represents  the  rate  of  sales  at  time 
The  lG3t  three  variables  tire  related  by  the  condition: 

(1.2)  H(t)  -  P(t) 

that  the  rate  of  increase  in  inventories  is  equal  to  the  rate  by  which 
production  exceeds  sales, 

Since  a  rationale  for  a  cost  formulation  of  this  sort  has  been  given  clna 
where  [2],  further  explanation  is  probably  not  warranted  here.  Limited 
experience,  however,  indicates  that  estimation  of  the  parameters  is  feasible 
and  that  the  decision  rules  derived  on  the  basis  of  such  estimates  are  only 
moderately  sensitive  to  errors. 

In  Section  2  we  xfili  be  concerned  with  finding  a  rule  for  setting  the 
production  rate  optimally,  in  the  sense  that  the  expected  value  of  the  cos~ 
function  (1.1)  is  minimised,  taking  into  account  new  information  available 
to  the  decision*- '•fiakev  which  is,  in  turn,  reflected  in  revisions  of  forecast^, 
of  future  sales.  It  will  be  shown  that  the  optimal  rule  for  this  relative.-!* 
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simple  cost  functional  is  as  follows  i 

(1.3)  P(t)  -  V^00  o'*0  s(t,o)  &  - 

where  S(t,0)  is  the  forecast  of  sales  made  at  time  t  for  time  t+O.  The 
production  rat©  scheduled  depends  first  upon  a  ‘present  value1  of  expected 
future  sales  as  well  as  the  initial  state  of  inventories  of  finished  goods. 
Through  the  latter  term,  corrections  of  errors  in  previous  forecasts  aro 
gradually  made,  Subsoquent  sections  will  be  concerned  first  with  the 
response  characteristics  of  the  rulo  (1.3)  for  three  simple  kinds  of  fore¬ 
casts.  It  will  bo  shown  that  the  rule  is  a  ‘low-pass’  filter  (cf. Simon  [6]) 
ar.d  that  inventory  adjustments  might  lag  behind  sales  changes  for  forecasts 
tending  to  extrapolate  present  sales  conditions  into  the  future.  The  sensi¬ 
tivity  of  the  rules  and  costs  associated  with  their  performance  to  errors 
in  estimating  the  coefficient  t  ,  in  the  3ales  forecasts,  and  in  forecast- 
ing  over  only  a  limited  ’horizon’  will  then  bo  examined. 

Although  the  work  steas  from  production  scheduling  problems  in  a 
specific  fine,  formulations  of  this  sort  appear  to  be  applicable,  with 
modifications,  to  other  situations.  It  is  possible,  too,  that  explicit 
statements  of  a  criterion  as  in  (l.l)  will  provide  some  clues  to  better 
design  of  physical  servomechanisms  (e.g. ,  control  of  temperature,  pressure, 
and  recycle  rates  in  a  cracking  tower  together  with  the  design  implications 
of  better  control).  Furthermore,  such  techniques  are  potentially  useful  in 
the  study  of  business  fluctuations-  In  particular,  the  derivation  of 
behavior  hypotheses  from  coet  criteria  not  only  provides  a  bridge  between 
the  theories  of  the  fira  and  business  cyclos,  but  also  suggests  the  structurs 
of  such  relations  for  statistical  estimation  (decision  variables,  data 
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relevant  to  these  decisions,  tho  role  of  expectations)  ar*i  provides  an 
independent  source  of  estimates  of  the  parameters  of  those  relations, 

2,  Derivation  of  Optimal  Roles 

Ths  problem  posed  before  may  be  restated  in  a  simpler  fora  by  extra¬ 
polating  the  * certainty -equivalence*  proof  of  Simon  [7)  to  the  situation 
in  which  decisions  may  change  continuously  in  time.  Let* 

P(t,9)  be  the  rate  of  production  planned  at  time  t  for  tini9  t<9, 

H(t,0)  be  the  level  of  inventories  planned  at  time  t  for  time  t+9,  ai-d 
S(t,Q)  be  the  rate  of  sales  forecasted  at  time  t  for  time  t+d. 

The  plan  which  minimizes  the  following  coot  functional  is  the  optimal  rule 
for  the  dynamic  problem: 

(2.?.)  C«15jii  )2  *  (P-P ')*]  dO, 

■S-a,  ?J0  •  *  ' 

where  ?  is  a  positive  number  and,  for  simplicity,  Hg  and  Pc  are  constants 
Wo  then  wish  to  find  a  way  of  minimising  (2.1)  subject  to  tho  following 
consistency  relation: 

(2,2)  H*  ■=  F  -  S, 

the  prime  indicating  differentiation  with  respect  to  9, 

The  problem  is  to  find  some  derivative  P^(t,0)  which,  when  set  equal, 
to  P^°^(t)  the  actual  production  rate  at  time  t,  will  minimize  the  expeotra 
valuo  of  equation  (l.l).  The  calculus  of  variations  will  yield  tho  conditions 
for  a  minivans  coot  in  the  planning  3pace.  With  tho  aid  of  the  Laplace  trrns 
fora,  we  will  then  solve  the  resulting  conditions  for  0  -  Theso  cqu&t 
represent  the  optisal  decision  rule.  Having  gone  through  those  icandpulaf-.-.-ti 
we  will  examine  cone  of  the  qualitative  pi-operties  of  such  rules,  perhaps 
gaining  insist  as  to  possible  applications  and  limitations  of  this  appro-rcl:-. 
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2.1  Conditions  for  Hininma  Coato  If  wo  lot  P*(0)  be  the  optimal 
production  plan  over  the  future  (wo  ore  dropping  the  argument  t  for  the 
time  being),  an7  other  plan  may  be  v/ritten  as  a  'variation'  from  this 
optimums 

(2.3)  P  -  P*  ♦  v». 

Kxcept  that  it  be  differentiable  and  satisfy  certain  end-point  conditions 
to  bo  specified  below,  the  function  v(0)is  arbitrary.  Sin^e 
H*  ■  (P*  -  S)  ♦  (P  -  P*)  -  H*  ♦  v', 

It  follows  that 

(2.4)  H  -  R»  ♦  v 

plus  a  constant  which,  as  we  shall  see  later,  is  equal  to  zero. 

Substituting  (2*3)  and  (2*4.)  into  the  functional  (2.1),  an  equation 
of  the  following  fom  is  oDtaincd: 

(2.5)  C  -  C*  ♦  2L  ♦  Q. 

the  avorage  cost  ueing  tho  policy  P*  ia  given  by: 

T 

(2.6)  C*  •*  11a  \f  t^2(H*-Hc)2  ♦  (P*-Pc)2  ]  dS. 

T-*co“  0 

L  represents  the  line.tr  conns  in  v  and  its  derivative: 

(2.7)  L  -  lio  \  f  [  32(H*~H  )v  +  (P*-P  )v']  dO 

T-w  TJQ  c  c 

and  Q  represents  the  quadratic  terms: 

(2.8)  Q  -  lim“/’T  C^W)2  ♦  (t*)2]  d©. 

T-coTy  0 

Sinco  P*  is  tho  minimum-cost  plan  for  production,  for  all  the 
admissible  functions  (alternatives)  v,  we  must  have: 

(2*9)  C  -  C*  ♦  2L  ♦  Q  £  C*, 
and  hence  that 
(2.10a)  L  -  0,  and 

(2.10b)  Q  >  0 
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for  &11  admissible  functions  v.  Because  the  ratio  of  cost  coefficients, 
o  ,  is  positive,  we  can  see  immediately  from  (2„8)  that  condition  (2.10b) 
for  a  minimum  io  satisfied* 

What ,  then,  are  the  conditions  under  which  L  ■  0?  Clearly,  this 
holds  for  the  trivial  case  ■  P*-Pc"  0,  which  implies  a  very  special 

kind  of  sales  forecast  (without  error).  P*  and  H*  may  be  found  under  more 
general  conditions  as  followa->  Equation  (2.7)  may  be  written  in  terms  of  v 
(and  not  its  derivative)  in  the  interval  (0,T)  and  ite  values  at  the  ond~ 
points  as  follows! 

(2.11)  L  -  lim  [#(H*-H  )  -  (P*')]v  d©  ♦  (P*-P  )v|T?-  0. 

The  intejral  will  vanish  if  the  following  differential  equation  in 

planned  future  production  rate3  and  inventcy  level?.  1j  auxirfied: 

(2.12)  tf2  (lf»  -  Kg)  -  P*‘  -  0. 

For  the  equation  to  be  satisfied,  it  is  also  necessary  to  specify  that 
v(0)  ■  v(T)  ■  0.  The  condition  on  the  ’variation'  implies 
(2.13a)  H*(0)  -  H(0), 

(2.13b)  H*(T)  -  H(T). 

That  is,  the  initial  and  terminal  inventory  levels  must  be  specified.  The 
first  condition  means  that  H(0)  is  part  of  the  data  upon  which  the  rate  of 
production  P*(0)  is  set.  We  shall  see  below  that  if  this  condition  were 
not  imposed,  it  would  be  possible  to  have  infinite  average  co3t3  per  unit 
time.  Tho  second  condition,  (2.13b),  will  not  bo  taken  seriously  because 
in  taking  the  limit  as  T  -*  oo ,  the  terminal  inventory  can  be  specified 
arbitrarily  without  affecting  tho  decision  for  the  immediate  future. 

2.2  Solution  for  Ivisediate  Decision.  Of  course,  tho  differential 
equation  (2.12)  is  the  classical  Euler-L a grange  condition  for  an  extremum 
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of  a  functional.  He  are  not,  however,  30  such  interested  in  finding  the 
time-path  P*(0)  as  in  the  initial  decision  P*(0).  Subsequent  plans  will, 
of  course,  be  made  on  the  basis  later  information.  To  9olve  the  differ¬ 
ential  equation  for  P*(0)  and  later  planned  production  rates  we  will  use 
an  integral  transformation  that  has  been  very  useful  in  the  analysis  of 
other  linear  system*. 

Th«  Laplace  transform  of  a  function  f  (©)  over  the  interval  (0,T)  is 
defined  as* 

✓T 

(2.H)  f„  (p)  -[  e**13®  f(»)  d». 

f  '0 

The  transformed  function  no  longer  depends,  of  course,  upon  9,  but  on  tho 
parameter  p,  a3  well  as  the  upper  limit  of  integration,  T.  fr(p)  may  be 
interpreted  as  tho  present-value  of  a  timo— series  f (0) ;  it  is  a  function 
of  the  interest  rate,  p,  and  not  time.  It  follows  from  integration  by 
parts  that* 

(2.15)  4k)  (p)  -  Pk  fT(p)  -  ^  p^Cf^^CO)  -  e~pT  f(j_l)  (T)  3, 
befinlng 

(2.16)  f (p)  -  lia^^  ^(p). 
it  follows  that* 

(2.17)  "?kT  (p)  -  Pk  f(p)  -  X  Pk_J  f(j"X)(0) 

as  long  as  f  (©)  and  its  derivatives  do  not  approach  infinity  core  quickly 
than  a  polynomial. 

Denoting  tho  transform  cf  P*(0)  over  the  interval  (0,T)  by 
etc. ,  we  will  tronoform  the  differential  equation  in  0  into  an  algobraic 
equation  in  p,  obtaining* 

(2.18)  (l-e~PT)3  -  p?T*  LP(0)  -  o~PT  P(T)]  -  0. 

Prom  the  definitional  relation  (2.2)  we  also  have: 

(2.18)  p  Kj,  -tH(0)  -  e'pT  H(T)]  -  PT  -  Sr 
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Substituting  (2.19)  into  (2.18)  and  rearranging  terms,  we  obtain: 

(2.20)  UV"1  -  p)^-  tfW**  [H(T)-J^]  -  e*^P»(T)  ♦  P»(0) 

-  -  tH(0)-flc)  . 

3  “  sj  exists  for  all  p  having  positive  real  parts  if  the  sales 

forecast  S(0)  is  bounded.  In  the  Unit  as  T  -*  co ,  the  terminal  inventory 
level,  H(T),  becomes  irrelevant  and  P*(T)  drops  out,  leaving: 

(2.21)  (ftV-1  -p)P  ♦  P*(0)  -  $2p~l  {s  -  [H(0)-Hc]J. 

The  inventory  level  at  the  beginning  of  the  period,  H(0),  and  the 

forecast  of  sales,  P(Q),  are  pieces  of  information  used  in  setting  the 
production  rate,  P*(0),  for  the  instant  of  time  immediately  ahead.  That 
H(o)  is  not  a  decision  variable  cooes  from  the  end-point  conditions  froa 
the  variational  analysis,  viz.  H*(G)  *  H(0)r  However,  P*(0)  is  not  so 
specified.  It  may  be  chosen  no  that  the  differential  equation  (2.12),  or 
equivalently  the  transform  equation  (2.21),  will  be  satisfied. 

P*(0)  may  be  found  as  follows.  The  'present  value*  of  planned  future 
production  rates,  P,  is  finite  for  all  p  having  positive  real  parts  aa 
long  as  P(O)  does  not  approach  infinity  any  faster  than  a  polynomial 
(ordinarily  it  would  be  bounded).  The  coefficient  of  P  will  vanish  for 
values  of  p  (say,  pjj)  satisfying  the  equation: 

(2.22)  i2  p"1  -  Pk  -  0 

as  well  as  the  condition  Re(p)  >0.  Consequently,  the  desired  root  is 
where: 

(2.23)  -  ♦  *  .  k»l,2,  respectively. 

Substituting  into  equation  (2.21),  we  then  obtain  the  expression 
for  the  optimal  production  rate  in  the  'next  instant*  of  time  as  a 
function  of  initial  inventory  levels  and  a  forecast  of  future  sales: 

(2.24)  P*(0)  -  i  S(o)  -  ^H(0)-Ho). 


t 


'£ 
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R ©-Introducing  the  time  argument,  t,  equation  (2.24)  nay  be  written  as* 

'00  -$) 

'0 

The  rule  states  that  the  production  rate  to  be  set  for  the  t*th  instant 


(2.25)  P(t)  -  P*(tB-»0)  -  e^  S(t,0)  dO  -  tftH(t)-H0). 


of  time  should  be  a  weighted  sum  of  future  expected  sales,  less  some 
partial  correction  of  inventory  levels  (arising,  in  this  case,  from 
previous  sales  forecast  errors).  We  can  see  from  the  equation  a  rationale 
for  the  stipulation  H*(t8«0)  “  H(t).  There  has  been  no  restriction  upon 

the  forecast  function.  Consequently,  if  H(t)  were  directly  set  by  the 
rule  and  there  were  a  discontinuous  change  in  the  forecast,  P(t)  “  A(t)  ♦ 
S(t)  would  become  infinite,  and  so  would  the  costs. 

2.3  Solution  for  Distant  Plans.  It  is  also  possible  to  calculate 
anticipated  or  planned  rates  of  production  In  the  future.  We  will  start  with 
equation  (2.20),  interpreting  T  not  a3  the  end  of  the  planning  horizon  as 
before  but  as  the  running  argument  in  the  'planning  space*.  Since  H*(t)  is 
no  longer  the  specified  end-point,  the  asterisk  is  used*  By  choosing 
as  ♦  (,  respectively,  the  identity  (2.20)  gives  the  following  two  equations 
in  H*(T)  and  P*(T) : 

(2.26a)  -*tH»(T)-Hj  -  P»(T)  -  6^  [ST(J0  -  S(tf)] 

(2.26b)  ftH*(l)*Hc)  -  P*(T)  -  -  e^lST(^)<*S(l()J  ♦  ZSs“^Eh(0)-Ho]. 


Eliminating  the  term  involving  inventories  from  this  system,  the 
following  expression  1 
(2.27)  P*(T)  -  | 

-df 


If*  4£T  ,S(0-T) 

iff  v  *0 

b  e-rf(0«T)  "| 

tH(0)-H  3* 

J 

dO 


1 


-  ID  - 


C 


For  large  values  of  T,  equation  (2.27)  may  be  approximated  bgrt 

(2  26)  P*(T)  -  rr  •  IS(T«0)  ♦  S(T-D)]  <». 

The  Inventory  adjustment  term  has  dropped  out  essentially  because  all 
the  adjustment  will  be  expected  to  take  place  over  a  relatively  short 
period  of  time,  Tince  the  calculation  Is  completely  in  the  ‘ expectations 
space'  there  will  be  no  new  'shocks-  from  the  forecast  errors  to  correct. 

2.4  Relation  to  the  Rules  Discrete  In  Time.  It  is  worth-while 
to  examine  at  this  point  the  relation  of  the  ‘rules'  presented  here  with 
the  work  of  Holt,  Modigliani,  and  Simon  12].  If  instead  of  varying 
continuously  in  time,  the  production  rate  were  set  for  a  period  of  length 
h,  the  cost  function  would  take  the  form  of  a  sum: 

(2.29)  C  -  lia  ?  L  1^%-H.)2  ♦  tt*2  (P^-P-)2] 

where  Hq  is  the  inventory  level  at  the  end  of  the  O'th  period  of  time  in 

the  future.  Pq  is  the  amount  of  production  during  the  O’th  period  of  length  h, 

and  the  summation  is  over  0  -  h,  2h,  ....  X  h. 

* 

The  decision  rule  which  leads  to  minimizing  the  expected  value  of 
the  cost  function  above  may  be  shown  tc*  be: 

(2-30)  Vh  -  j”  '•i" 3t.s  - 

*1 

The  parameter  ^  is  related  to  the  oost  ratio  as  follows: 

(2.31)  ^  -  \  [2  ♦  X2*2  - 

If  the  period  is  rather  short,  the  parameter  Xj  is  approximately 
1  -  hjf.  Consequently  the  production  rule  may  be  written  as 
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(2.32)  pt^,  -  W  a  -  bi0<O'h)/h  \.e  -  **  (W 

As  wo  lot  h  approach  zero,  wo  have 

1  oo  (Q-h)/h  l  . 

(2.33)  p(t)  ■;  Pw  ■?  t  (1  -  h»  g  St>0  h  -  S  (Hf-H,) 


*fo 


Oh 

CD  -fa  ,  . 

*  3(b,fr)  Od  -  6  J» 


10  ~  ~  o' 

This  Is,  of  course,  identical  with  the  continuous  rule  (2a25)  which  has 
been  previously  obtained^ 

If,  on  the  other  extreme c  we  let  the  time  interval  become  very  large p 

-2 

then  is  approximately  equal  to  (ho)  .  Hore  the  rule  would  become 

(2.24)  pt<h  -  u-h-V*  )  X,  <*>#>'2Ce*1)  -  &*>»{■•  j 

~  St.h  -  -  »o>- 

This  limit  of  tho  production  adjustment  which  is  due  to  Metzler  (A]  and 
analogous  to  that  of  Vassian  l8j,  can,  of  course,  be  obtained  from  a  cost 
function  involving  inventory  costa  alone. 

3«  Dynamic  Performance  Characteristics 

3c  1  F.esy.ae  to  a  Sinusoidal  Input 0  Returning  now  to  the  deoision 
rule,  aquation  (2.25),  it  is  evident  that  the  dynamic  performance  will  depand 
critically  upon  ths  type  of  forecast  that  is  used.  We  will  examine  the 
performance  of  the  rule  for  a  sinusoidal  sales  pattern  for  three  types  of 
forecasts! 

1.  'Perfect':  S(t,0)  -  S(t+0), 

2.  'Null':  S(t,C)  -  0, 

3.  'Naive'*  S(t,0)  -  S(t). 

Other  possibilities  involve  mistakes  as  to  amplitude,  frequency,  and 
phase  of  the  fluctuations,  but  we  will  not  consider  these  here. 


With  the  porfeot  forecast,  the  production  performance  may  be  obtained 


directly  froa  the  Euler-Lagrange  equation  (2.12)  to  box 
(3.1)  S2  P(t)  -  P(t)  -  ^(t). 

The  relative  magnitude  of  production  fluctuations  to  sales  fluctuations 
is: 

(3*2)  lrrl  -  rr -fer 

where  o>  is  the  angular  frequency  of  fluctuations  of  sales.  The  lag  in 
production  adjustments  behind  eales,  ftp,  would  be  zero,  as  we  would 
expect  with  a  perfect  forecast.  The  relative  response  jYpjie  close  to 
unity  for  &><}{,  but  drops  off  eharply  for  large  a).  In  other  words, 
production  adjusts  to  slow  changes  in  sales,  but  ignores  those  of  a  very 
short  duration.  The  'cut-off1  frequency  depends,  of  course,  upon  the 
relative  costs  of  Inventory  storage  and  overtime,  which  are  reflected  In 
the  coefficient  The  relative  magnitude  of  inventory  adjustments  to 
the  fluctuations  in  sales  is  given  by: 

(3.3)  M 

and  /Zp  ■  tt/2  (i-o.  the  inventory  reaches  its  mnrlmraa  level  when  production 
oquals  sales  and  both  are  increasing  The  function  }Zp|is  small  for  very  hi^i 
and  very  low  frequencies,  reaching  a  marl  roan  of  1/2&  for  X» 

The  analysis  with  the  perfect  forecast  would  give  tho  beat  pcrforaar  oe , 

It  ie,  however,  dynamically  unstable,  so  that  any  errors  in  the  fore  caste 
would  ultimately  'blow  up*.  It  can  be  shown  that  the  performance  of  the 
rule  is,  in  general,  dynamically  stable  if  the  forecast  ie  what  servo¬ 
mechanism  engineers  call  'physically  realizable'  (i.e.  depends  only  upon 
past  data). 
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Perhaps  the  most  nipple  realisable  forecast  is  the  'null*  forecast, 
that  future  sales  are  expected  to  be  equal  to  some  long-term  are  rage. 

With  S(t,0)  ■  0,  the  rule  (2.25)  beooaee* 

(3.4)  P(t)  -  H(t)-H0). 

Consequently  the  relative  magnitude  of  production  and  sales  fluctuations 
is* 

<3-5>  N 

Since  all  values  of  the  function  are  less  than  unity,  the  rule  with  the 
null  forecast  responds  more  to  sales  fluctuations  than  that  with  the 
perfect  forecast  for  all  frequencies.  The  angle  by  which  production  lags 
behind  sales  is  given  by* 

(3.6)  Aq  -  tan  0, 

increasing  with  the  frequency.  The  magnitude  of  response  of  the  inventory 
level  ie: 

(3.7)  |ib|*  -i  -—==• 
with  the  same  lag  as  production, 

A  'naive*  forecast  is  equivalent  to  the  stipulation  that  S(t,0)  -  S(t) 
This  leads  to  the  production  rule: 

(3.8)  P(t)  -  S(t)  -tf[H(t)-Hc). 

In  a  way,  this  case  is  rather  uninteresting,  because  production  will  equal 
sales  and  consequently  there  will  be  no  fluctuations  in  the  level  of 
inventories, 

3.2  Sales  Impulse  Anticipated  T  Units  of  Time  In  Advance,  Another 
way  of  examining  the  offects  of  the  rule  is  the  following.  Suppose  sales 
and  production  have  been  maintained  for  a  long  time  at  a  constant  rate 
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(which  can  be  assumed  equal  to  zero  without  loss  of  generality).  At  time 
t  "  T,  an  impulse  of  sales,  which  was  anticipated  at  time  t  ■  0,  occurs. 
It  is  then  possible  to  examine  the  cost  of  production  as  a  function  of  Tt 
the  length  of  time  in  advance  that  such  a  sale  is  expected. 

The  sales  function  would  then  be 


(3.9)  S(t)  -  6  (t-T) 
where  6  is  the  l 

(3.10)  S(t,0)  - 


where  6  is  the  Oiruc  delta.  The  sales  forecast  function  would  be 

{b(t*©-T),  0  <t<T,  0<0, 

0  ,  otherwise. 

If  T  »  0,  there  would  be  no  advance  notice  of  the  sale)  if  T  ■  co. 


there  would  be  perfect  advance  knowledge.  These  two  extremes  have  been 
examined  by  Holt  and  Simon  [3]  with  -different  methods  than  those  to  be 
used  here. 

With  some  reasonably  convenient  operations  involving  the  Laplace  trans¬ 


form.  equations  (3.9)-(3.10)  and  the  decision  rule  lead  to  the  following 
function  for  the  production  'response' ; 

ro  *  t  <  o, 

(3.11)  P(t)  -  ]  Ye“*T  cosh  C  ',  t  <  T, 

(^e“at  cost  tfT-  ▼  <  t.~ 

Similarly,  the  time-path  of  the  Inventory  levels  may  be  found  by  integration 


to  he: 

(\ 

t 

t  <  0, 

(3.12) 

H(t)  -  ).’** 

1-tT ^ 

sinh  ^t, 

sinh 

0  <  t  <  T 

T  <  t. 

Substitution  of  (3.11)  and  (3.12)  into  the  cost  function  (l.l)  leads 
to  the  total  cost: 

(3.13)  C(T)  -  X  (1  ♦  *”2^T)/2. 


o 
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Evidently  C  is  a  decreasing  function  of  T  and 
C(0)  -  t ,  C(oo)  -  Y/Z, 

The  'success'  of  the  forecast  is  defined  ast 


(3.14) 


This  function  measures  the  fraction  of  total  cost  savings  In  forecasting 
can  be  realized  with  a  forecast  having  horizon  T.  Fran  the  cost  function 
(3.13).  the  'success'  function  becomes: 

(3.15)  s(T)  -  1  -  s’2*  . 

With  •  .09*  a  forecast  horizon  of  5  months  would  be  9556  'uuccessful.' 
Although  the  entire  future  is  relevant  to  some  extent,  the  cost  effect  of 
errors  declines  for  this  case  quite  rapidly.  If  the  co'-t  of  inventory 
deviations  is  higher,  only  shorter  horizons  would  be  required;  as  the 
storage  cost  completely  dominates  the  cost  of  deviations  in  the  production 
rate,  only  the  Immediate  future  is  relevant  (as  in  the  Hetzler  model  [4]) 

4.  Effects  of  Forecast  and  Estimation  Errors 

Since  Imperfections  in  the  information  available  to  the  decision¬ 
maker  is  an  important  characteristic  of  planning  problems,  it  is  certainly 
relevant  to  examine  the  effects  of  various  kinds  of  errors  upon  not  only 
the  decisions  that  are  made  but  also  the  extra  costs  that  are  incurred  as 
a  result  of  these  errors.  Possible  sources  of  error  in  the  types  of  models 
considered  here  might  be  summarized  as  follows: 

1.  Estimation  of  parameters : 

i.  Cost  coefficients  (**).  .  . 

ii.  Forecasts  (S(t,0)  -  S(t*0)  ), 

ill.  Feedback  information  (H(t,0)  -  H(t))j 


© 
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2.  Performance  (e.g. ,  differences  between  planned  production  and 
actual  production,  P(t,0)  -  P(t))| 

3*  Specification  of  the  criterion  function: 

I.  Tine  grid:  planning  horizon  (T)  and  decision  frequency 
(l/h  for  the  rule  discrete  in  tlno,  infinite  for  the 
continuous  rule) , 

II.  Aggregation  of  components:  product  of  the  firn  an.  the 
organizational  units  and  locations, 

ill.  Excluded  components:  decision  variables  (e«g« (  raw  materials 
purchases,  vork-in-proceas  Inventories,  etc.)  and  coats, 
iv.  More  fundamental  assumptions:  the  quadratic  model,  lack  of 
unllatorial  restrictions,  the  expected  value  of  the  function 
as  a  criterion  for  choice,  etc. 

Although  little  can  yet  be  said  about  the  sources  of  error  toward  the 
bottom  of  the  list,  some  statements  can  be  made  about  the  effects  of 
cost  estimation  errors  and  Information,  control,  and  forecast  errors* 

This  we  do  below* 

&.1  Errors  In  Cost  Coefficients.  Since  the  parameters  of  the  decision 

rules  depend  upon  the  values  of  the  cost  coefficients,  the  behavior  indicated 

will  be  costly  to  the  extent  that  incorrect  estimates  of  the  coefficients 

were  used*  In  the  case  involving  inventory  and  production  rate  ooots,  an 

error  in  the  ratio  ^  and  the  parameters  H  and  P  will  result  in  different 

c  o 

* 

decisions  (and  hence  higher  costs)* 


*  Costing  errors  in  decision  rules  was  first  carried  out  by  Roberson, 
Holt,  and  Modigliani  [5]  for  the  rule  based  on  inventory,  overtime,  and 
labor  turnover  costa0  H.  Theil  has  made  extensive  investigations  for  a 
static  quadratic  welfare  function  in  work  that  is  as  yet  unpublished* 
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The  cost  functional  (1*1)  nay  be  written  ass 

(4.1)  c-  i2  TiwTJ2  ♦  (p-p  )2< 

O  6 

The  bare  Indicate  the  tine  average  of  the  square  of  Inventory  end  production 
deviations *  respectively* 

Assume  that  sales  can  be  forecasted  perfectly*  that  there  is  no  error 
in  estimating  Hg  (Pc  Is  irrelevant)*  and  that  sales  may  be  represented  by 
the  undamped  sinu soldi 

(4.2)  S(t)  -  S  ♦  A  cos  a  t,  A  •  const. 

The  resulting  production  rates  (If  sales  are  perfectly  known)  will  obey 
the  following  differential  equation! 

(4.3)  ^P(t)  ♦  o2  P(t)  -  c2  S(t) 

If  c2  is  the  estimate  of  tf2  .  The  production  performance  would  therefore  be 
given  by* 

(4.4)  P(t)  •  A  oos  tit  ♦  S 

—  i.e. ,  production  would  have  the  sane  frequency  and  phase  as  sales  and 
an  amplitude  no  greater  than  that  of  sales  whatever  the  value  of  the 
estimate  of  the  cost  ratio. 

The  time  average  of  production  fluctuations  would  then  be! 

(4.5)  (M,)J  -  4 

Similarly,  it  can  be  shown  that  the  average  of  the  inventory  fluctuations  is: 

(4.6)  (Ml.)**  L  • 

If  c2  were  too  large*  then  more  of  the  fluctuations  In  sales  would  be  ’taken 
up’  in  production  fluctuations  than  would  be  deisrable  (and  conversely)*  but 


neither  of  these  errors  in  the  estimate  of  the  cost  coefficient  Mould 
sake  much  difference  at  low  frequencies  (i.e.  if  the  sales  rate  changes 
slowly). 

Assuming  that  S  -  P„,  the  time  average  of  costa  can  bo  written  as  a 


function  of  the  true  coefficient  its  estimate  <? ,  and  the  frequency  of 
sales  fluctuations  from  (4»l)  as  follows! 


Similar  results  would,  of  course,  bo  obtained  with  a  rule  with  which 
decisions  are  made  at  distinct  points  of  time.  The  relative  Increase 
in  costs  due  to  an  error  of  estimating  the  cost  coefficient  1st 


Coln  tfj 

If,  for  example,  the  coat  ratio  c^  were  3.056  greater  than  ■  .3 

per  month  (as  is  approximately  the  case  for  the  factory  studied  in  (2)) ,  and 


w  ■  „525  radian s / month  (one  cycle  per  year),  then  the  error  would  increase 
costs  by  &  factor  of  about  0,1$6.  It  is  a  distinct  advantage  that  the 
performance  of  the  rule  and  the  coat  of  the  performance  is  not  very 
sensitive  to  errors  in  the  estimates  of  the  coefficients,  (Extreme 
insensitivity  is,  however,  equivalent  to  stating  that  costs  are  not  affected 
much  by  the  scheduling  decisions.)  Furthermore *  substantial  deviations  from 
predictions  of  inventory  fluctuations  made  with  rules  of  this  sort  might  be 
expeetod  because  coot  differences  might  be  hard  for  firms  to  perceive. 


Equation  (4.6)  is  obviously  a  function  of  the  frequency  of  the  sales 
fluctuations.  One  can  see  that  the  costs  ’peak'  at  the  angular  frequency 
<a  m  c,  the  square-root  of  the  estimated  cost.  The  maximum  value  of  the 


o 
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relative  cost  difference  is: 

(^.9)  i  (t-f)  . 

Thus,  the  maximum  relative  increase  in  cost  for  c2  -  1.10  would  be 

0.23$. 

4.2  Forecast,  Control,  and  Information  Errors.  Letting  h  •  1,  the 
discrete  form  of  the  rule,  (2.30),  becomes: 


(4.10) 


-V» 


whore  Pt,l  *8  amount  of  production  planned  for  the  coadng  period,  H^  q 
if  the  (imperfect)  estimate  of  inventories  at  the  end  of  the  proceeding 
period,  and  ^  is  the  forecast  made  at  time  t  for  sales  in  the  (t-*0)lth 
period.  If  there  were  no  errors  in  these  variables,  they  would  be  equal 
to  P^4^,  „  and  S^^q,  respectively. 

Suppose  that  the  time— series  is  stationary  and  independent: 

-  -  — 2 

(4.11)  St  -  S  ♦  ,  *5  -  0,  *3  >  0, 


and  that  the  forecast  of  future  sales  is: 

(4»12)  St  q  ■  S. 

We  will  further  stipulate  the  'information*  error: 

(4.13)  Ht  0  -  Ht  ♦  »,,t  .  ^  *  0.  *h  >  °* 

and  the  ’control'  error:  _ 

2 

(4*14)  P^^  "  4  »  fcp  “  0*  *p  >0* 

We  will  assume  that  the  errors  are  statistically  independent. 

Actual  production  in  the  (t+l)'ot  period  is  than  described  by: 

(4.15)  Pt.x  -  -  (1-^)0^  .  -  h0)  .  s. 


O 
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The  difference  equation  above  can  readily  be  solved  with  the 
aid  of  the  power-series  transforms 


(4»  16) 


ffl  '  (1"V  Jc  ‘l'p.t-o  -£p.w  *  (1-S.>  £>*?  5S.W) 

(l-X, )  CD  ©  _  ’ 

*  “x^  xleV,t-o*<1’xi)tlH  4  3- 


From  the  relation  Ht  -  Ht-i  ■  Pt  -  St*  '!SO  can  also  finds 

(<.•17)  Ht  -  £  £  tp>W)  -  (M^  £  (V|, 


Consequently,  the  moan— square  deviation  of  production  is: 

(4.16)  (P-P0)2  -  __2_  rg  «.  2(1-Xx)2  X  l^^T  ♦  (3  -  P  )2. 

l^L  P  1  ♦  Xx  H  1+XX  ~ 

The  deviation  in  inveutorlea  is: 


(4.19) 


«*V 


"77 


4  am  *  3  • 


1-xJ  >* 


The  average  cost  of  thsae  errors  will  then  be: 
_  (l-x,'2  -  - 

(4.20) 


C  -  — (H-Hq)2  ♦  (p-pc)2 

.  1  ^  »  (1*X1)  ^  ♦  i^l  ^  ♦  (3  -  Pc)  . 

h  ?  *x  11  \  3 


Since  X^  ie  always  leas  than  unity  (and real),  it  ia  possible  to 
draw  certain  conclusions  about  the  effect  of  the  three  errors  upon  average 
costs.  The  rule  can  clearly  do  reasonably  well  with  sloppy  information 
about  existing  inventory  levels,  while  deviations  ftaei  planned  production 
are  rather  costly.  Although  the  model  is  not  particularly  sensitive  to 
forecast  errors  if  they  are  serially  independent,  the  rule  is  quite 
sensitive  to  a  bias  In  the  forecast. 


0 
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For  example,  if  o  ■  .09*  as  In  the  previous  illustrations ,  would 

be  equal  to  .74*  and  the  cost  function  (4.20)  would  become  i 

(4.21)  C  -  1.35  ♦  .09  ♦  .35  fcg  ♦  (S  -  P0)2. 

Of  course,  the  analysis  of  the  effects  of  error  is  by  no  Beans 
complete.  We  have  examined  their  effect  for  a  rather  simple  cost  function  — 
and  even  here  for  special  cases.  To  have  done  otherwise  would  have  led  us 
Into  too  such  analytical  underbrush.  However,  such  analyses  do  help  answer  the 
question  of  what  constitutes  a  good  forecast.  Obviously,  one  cannot  tell 
unless  he  has  some  idea  of  the  effects  of  errors  upon  the  criterion  for 
choice  (e.g.  production  costs).  The  costs  as  we  have  formulated  thm  are 
quite  insensitive  to  moderate  errors.  Again,  interpreting  the  ’decision 
rules*  as  behavioral  hypotheses,  the  insensitivity  of  costs  to  small  errors 
suggests  ^aviations  of  a  fair  magnitude  about  predicted  output.  If  in 
addition  a  'threshold  level*  in  cost  perceptions  is  postulated,  it  may, 
however,  be  possible  to  estimate  the  error  of  the  prediction,  at  least  at 
the  microeconomic  level. 
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SAID  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY  >'^  '7 
IMPLICATION  OR  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OR  ANY  OTHER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE, 
USE  OR  SELL  ANY  PATENTED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 


